Localized nonlinear waves in systems with time- and space- modulated nonlinearities 
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Using similarity transformations we construct explicit nontrivial solutions of nonlinear Schrodinger 
equations with potentials and nonlinearities depending on time and on the spatial coordinates. We 
present the general theory and use it to calculate explicitly non-trivial solutions such as periodic 
(breathers), resonant or quasiperiodically oscillating solitons. Some implications to the field of 
matter- waves are also discussed. 
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Introduction.- The concept of solitons was introduced 
in the seminal paper by Zabusky and Kruskal [H to char- 
acterize nonlinear solitary waves, which emerge from the 
delicate balance between nonlinearity and dispersion, and 
preserve their localized shapes and velocities during prop- 
agation and after collisions. It was soon realized that 
solitons are remarkably generic objects arising in differ- 
ent contexts and were observed in many experiments, for 
example in nonlinear optics starting with the pioneering 
works |2i]. Being solutions of integrable nonlinear wave 
equations, standard solitons appear usually in transla- 
tionally invariant systems involving nonlinearities which 
depend only on the relevant fields and do not have ex- 
plicit dependence on the space and time variables, be- 
cause such dependences typically break the symmetries 
of the system. However, those extra dependences of the 
nonlinear interactions open many new practical possibili- 
ties for generation and control of solitons, what has stim- 
ulated their study in the last thirty years. It was already 
recognized from the very first papers that specific de- 
pendencies of the equation coefficients on the spatial spa- 
tial coordinates and on time can preserve the integrabil- 
ity. Moreover this gave an origin to new concepts, like 
for example, non- autonomous solitons [^] highly relevant 
in nonlinear fiber optics. During the last few years these 
studies have seen an enormous revival in the context of 
systems ruled by nonlinear Schrodinger (NLS) equations 
because of their applications in the mean field theory of 
Bose-Einstein condensates (BECs) [HI, ^ and in nonlinear 
optics 0,0,13. 

In BEC applications, the matter waves being the nat- 
ural outcome of the mean-field description [5] have at- 
tracted a great deal of interest in the last years 0. 
The possibility of using Feschbach resonances to control 
the nonlinearities (see e.g. [o', To]) has lead to the pro- 
posal of many different nonlinear phenomena induced by 
the manipulation of the scattering length either in time 
111, ll2|, ll3| or in space [14]. In nonlinear optics, recent de- 



velopments [iG^ have led to the discovery of a new classes 
of waves such as the so-called optical similariton which 
arise when the interaction of nonlinearity, dispersion and 
gain in a high-power fibre amplifier causes the shape of 
an arbitrary input pulse to converge asymptotically to a 
pulse whose shape is self-similar. 

In this letter, we go beyond previous studies con- 
sidering space and time dependent nonlinearities and 
constructing different types of explicit solutions includ- 
ing: breathers, resonantly oscillating and quasiperiodic 
solitary waves. We provide explicit expressions for all 
those solutions, exploring realistic choices of nonlineari- 
ties which are experimentally feasible due to the flexible 
and precise control of the scattering length achievable 
in quasi-one-dimensional BECs [15] with tunable inter- 
actions. 

General Theory.- We consider physical systems ruled 
by the NLS equation 



(1) 



for an unknown complex function ?/^(t, x) with a potential 
v{x^t) and nonlinearity g{x^t). We focus on spatially 
localized solutions for which lim|a;|^oo ^) = 0- Out 
goal is to reduce Eq. ([T]) to the stationary NLS equation 



(2) 



where <I> = ^{X)^ X = X{t^x) is a function to be de- 
termined, /i is the eigenvalue of the nonlinear equation, 
and G is a constant. In this letter, we explore the case 
G = — 1, i.e., focusing or attractive nonlinearity. The 
case G = 1 does not pose any new challenges and will be 
studied elsewhere. 

To connect solutions of Eq. ([T]) with those of Eq. ([2]) 
we will use the transformation 



ilj{t,x) =p(t,x)e^^(''^^^(X(t,x)). 



(3) 



Other similarity transformations have been studied for 
NLS equations in different contexts 
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FIG. 1: [Color online] Solution (|17|) . for n = 1, cjo = 2, and 
e = 0.5. Initial data for Eq. {[11]) are i/(0) = \f2 and vt(fS) = 0. 

(a) Pseudocolor plot of |'0(x,t)|^ on the domain x G [—3,3] 

(b) Width v(i) (blue Une) and amplitude (red Une) vs t. 



Requiring ^{X) to satisfy Eq. ([2]) and ijjit^x) to be a 
solution of Eq. ([T|) with potential v{t^ x) and nonlinearity 
g{t^ x) we get the set of equations 

m+(pVx), = 0, (p'X,)^=0, X, + 2(p,X, = 0. 

Let us take 

X(t,x)=F(0, with =7(t)x + (5(t), (4) 

where ^{t) is a positive definite function representing 
the inverse of the width of the localized solution, and 
—S{t)/'y{t) is the position of the center of mass. After 
some algebra we obtain 

.A 
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p{t,x) 
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(p{t,x) 
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(5) 



(6) 



where a{t) is an arbitrary function of time. 

Thus, choosing 6{t)^'y{t) and F(^) (or equivalently 
p(x, t)) we can generate pairs g{x, t), v{x^ t) for which the 
solutions of Eq. ([T]) can be obtained from those of Eq. 
([2]) using Eqs. ([3|). We will use this fact to construct 
soliton solutions with an interesting nontrivial behavior. 

Localized nonlinearities.- As a first application of the 
method, we will study the case of localized nonlinearities. 
While the choice of possible models is very rich, because 
of applications the possible applications of our results to 
BECs with optically controled interactions [10.], we will 
focus on the gaussian shaped nonlinearity 



g{t,x) = -j{t)exp{-3e 



(7) 



and the potential including a combination of harmonic 
and dipole traps, which are typical in BEC experiments 

v{t,x) =cj^(t)x^+/(t)x + /i(t)-/i7^(t)exp(-2^^), (8) 



where 
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FIG. 2: [Color online] Width (blue line) and amplitude (red 
line) versus time of the solution of Eq. ([1]), given by Eq. (pT)) 
for n = 1, cjo = 2, and s = 0.5. Initial data for Eq. (|lip given 
by iy{0) = V2,iyt{0) = -6.125^2. 



Our choice corresponds to p{t,x) = y^exp(^^/2). 
In the particular case where J(t) = 0, /i = and a{t) = 
/ 7^(t)(it, we obtain 

v{t,x) =co\t)x^, g{t,x) = --i(t)e-^^'^'>\ (10) 

Thus, we have a harmonic potential and a modulated 
in time gaussian nonlinearity which would correspond in 
BEC applications to a modulated laser beam controlling 
the interactions optically. In this paper we focus on non- 
linearities and potentials given by Eq. (p!Q|) but many 
other possibilities can be explored using the same ideas. 
Defining = 1 /7 we rewrite Eq. ([9a|) in the form 



(11) 



which is the Ermakov-Pinney equation [19]. To discuss 
explicit solutions we choose (jo'^it) of the form 



a;^(t) = 1 + £ cos(a;ot). 



(12) 



where e G (—1,1) and 7^ cjq G M. After some algebra 
[19J , we obtain from Eqs. ^ and (fT2|) 



^(t)=[2yl{t)^2{y2(t)/Wf 



-1/2 



(13) 



where yi,2(i) are two linearly independent solutions of 
the Mathieu equation 



y„+4w2(t)y = 



(14) 



and W is their Wronskian. 

Thus it is Eq. ([T4j) that determines the temporal dy- 
namics of the obtained solutions. Considering different 
choices of the parameters cjq and £, one can single out 
three different types of behaviors which can be classified 
as follows: (i) periodic^ when £ = or in the frontiers 
bewteen the stability and instability regions of Eq. ([T4|) , 
(ii) resonant^ when £ 7^ and ^1,2 (^) belong to an insta- 
bility region of Eq. ([Hj), and (iii) quasi-periodic^ when 



3 




FIG. 3: [Color online] Plots of |'0(x,t)|^ for breathing solu- 
tions of Eq. ([1]), given by Eqs. {[TTl), for s = in Eq. ([12]), 
and 7(t) given by Eq. ([19}, corresponding to (a) n = 1 (b) 
n = 2. For both cases, x G [—6,6], t G [0, 10] and the initial 
data for Eq. ([TT]) are iy{0) = ^2, z/t (0) = 0. 




FIG. 4: [Color onhne] Solution of Eq. ([T]), given by Eq. ([TT]), 

for n = 1, and uj'^{t) = 1 + £cos(\/2t), with e = 0.5. The 
initial data for Eq. ([TT]) are z/(0) = 2^/^ and zyt(O) = 0. (a) 
Pseudocolor plot of |'0(x,t)|^ (b) Width zy(t). 



yi,2{t) are in the stability region. Below we give examples 
of all these cases. 

Concerning the solutions of Eq. ([2]), since /i = 0, Eq. 
([2|) becomes + = 0, and its solution is given by 

^X) =r]sn{r]X,h)/2dn{r]X,h), (15) 

with arbitrary r] and /c* = 1/ >/2. We choose X as 



/ 

J — ( 



(16) 



which takes values in (0,y^). This choice of X en- 
sures vanishing of ^ when x — oo. The imposi- 
tion of zero boundary condition for x +oo leads to 
T] = 2nK{k^)/^/7r where K{k) is the elliptic integral. 
This leads to the family of solutions 

1 



V^n = nV7(t)exp ( ^7^(^)^^ i(p{t, x))^i{nO{t, x) 



(17) 

where n = 1, 2, (p{t, x) = -(7^/47)^2 + / j'^{t)dt, 

^i{0) = y^K{h)sn{0,h)/dn{0,h) (18) 



and 0{t,x) = (2K{K)/^)p_^^U^e 
see that ilj{x^t) when \x\ ^ 00 in 
correspond to localized solutions. 

Resonant solitons.- The response of driven nonlinear 
systems to parametric perturbations is a field of interest 
which has recently also been investigated in the frame- 
work of the nonlinear dynamics of BECs [l^, 13, 21|. 

Starting with the particular case cjq = 2, in Eq. ([T2|) , 
we observe that for e 7^ we are in the instability region, 
thus, Eq. ([T7|) describes multisoliton resonant solutions^ 
since tpnix^t) possesses n — 1 zeros. In Fig. [H^a), we 
plot a resonant solution corresponding to n = 1, with 
ujQ = 2. Such resonant soliton has an increasing resonant 
behaviour. In Fig. (Tj^b), we plot its width and ampli- 
tude versus time. In Fig. [2l we also plot the width and 
amplitude versus time for the same soliton but now with 
transitory decreasing resonant behaviour, for n = 1. 



^ . It is easy to 
(fT7|) . thus they 



It is remarkable that we are able to construct explicit 
resonant solitons in a very complicated scenario including 
a localized pulsating nonlinearity and a parametrically 
modulated trap along the longitudinal direction. Our so- 
lutions provide the first explicit resonant solution in a 
parametrically modulated one dimensional BEC. A sim- 
ilar resonant behavior was observed in Ref. [21I] when 
driving a BEC, with positive scattering length, along the 
transverse direction. 

We have also studied the stability of these multisoliton 
solutions given by Eq. ([T7|) by computing numerically 
their evolution under finite amplitude perturbations. We 
have found that n = 1, which corresponds to the ground 
state, is a stable solution while higher solitons are un- 
stable (small perturbations take the solution far from its 
original profile). 

Breathing solitons.- When s = 0,a;(t) = 1, and then 



7(t) = V^[l + 3cos^(2t)] 



-1/2 



(19) 



is periodic. Then, the solutions given by Eq. ([T7|) are 
explicit multisoliton breathing solutions. In Fig. [3l we 
plot some of them corresponding to n = 1,2. Only for 
n = 1 (the ground state) the breathing soliton is stable. 

Quasiperiodic solutions.- In order to give an example of 
a quasiperiodic solution of Eq. ([1]) we choose uq = \f2 in 
Eq. ([T2|) ensuring that the solutions ^1,2 (^) of the Math- 
ieu equation ([T4j) belong to the stability region. Then 
Eq. ([T2|) has two incommensurable frequencies. In this 
way, 7(t), is a quasiperiodic solution and the solutions 
(fTTj) show a quasiperiodic behaviour. An example of this 
behavior is shown in Fig. (H Using Eqs. ([TTj) . we can 
also construct multisoliton quasiperiodic solutions. 

Moving solitons.- As a final example of the many possi- 
ble solutions which can be constructed using this method, 
we will present solutions of Eq. ([T]) when the centre of 
mass of the soliton moves with non-zero velocity. To do 
so, we set /(t), in Eq. ([9b|) equal to zero, with (5 7^ 
which leads to the following equation for S{t)^ 



5u-2{-fth) Jt+ 47^(^ = 0. 



(20a) 
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FIG. 5: [Color online] Pseudocolor plots of solutions of Eq. 
([T|), when the centre of mass of the soliton moves with non- 
zero velocity, (a) Quasiperiodic solution, with x G [—3,3], 
t e [0,40], for s = 0.6 and z/(0) = (2y2)^/^ zyt(O) = 0. (b) 
Resonant increasing solution, with x G [—6,6], t G [0,40] for 
e 0.3 and z/(0) ^2, z^t (0) 0. 



In this situation, the potential v and the nonlinearity g 
are still given by Eq. (p!Q|) by just taking 



at=7'(l+<5')-5?/47'- 
Fortunatelly, Eqs. ([2Q|) can be solved to get 



(20b) 



= 5^ cos(2r(t)), = r(t) + sin(4r(t)), (21) 

where r{t) = J j'^{t)dt^ and is an arbitrary constant. 
Hence, we can again obtain breathers, resonant solutions 
and quasiperiodic solutions while, at the same time the 
center of mass of the soliton moves in a complex way 
according to Eq. ([2T]) . As examples, in Fig. O we show 
quasiperiodic [Fig. [5^a)] and resonant solitons [Fig. [SJb)] 
while the center of mass moves according to Eq. (|2T]) . 

Conclusions.- We have constructed explicit solutions 
of nonlinear Schrodinger equations in complicated yet re- 
alistic scenarios with spatially inhomogeneous and time 
dependent potentials and nonlinearities. Although the 
range of nonlinearities and potentials for which this can 
be done is very wide, we have focused our study in mod- 
ulated harmonic potentials and localized nonlinearities. 
We have explicitly calculated breathers, resonant soli- 
tons, quasiperiodic solitons and solutions with nontrivial 
motion of the center of mass, related to different situ- 
ations of physical interest in the field of Bose-Einstein 
condensates. The ideas contained in this paper can also 
be applied to study multicomponent systems, nonlinear 
optical systems, higher-dimensional situations, etc., and 
may help in the design of potentials and nonlinearities to 
control the dynamics of Bose-Einstein condensates. 
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